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Functional impact of a recently identified quantitative 
trait locus for hippocampal volume with genome-wide 
support 

S Erk 1,2,12 , A Meyer-Lindenberg 3 ' 12 , P Schmierer 1,2 , 0 Grimm 3 , H Tost 3 , T Muhleisen 4 , M Mattheisen 5 , N Seiferth 1 , S Cichon 6 , 
M Rietschel 7 , MM Nothen 8 ' 9 , A Heinz 13 and H Walter 1 ' 2 ' 10 ' 11 ' 13 

In a large brain-imaging study, a multinational consortium has recently identified a common genetic variation in rs7294919 being 
associated with hippocampal volume. Here, we explored whether this quantitative trait locus also affects hippocampal function 
using a previously established reliable neuroimaging paradigm. We observed a significant effect of rs7294919 variation in the 
right hippocampus showing that hippocampal activation increased with the number of risk alleles. Furthermore, the risk allele 
was associated with decreased performance in a verbal learning and memory task. By showing that this single-nucleotide 
polymorphism also relates to behavioral difference and underlying brain activation in memory, our findings support the idea that 
rs7294919 may affect the individual capacity to resist disease in terms of diminishing or boosting hippocampal resources. 
Translational Psychiatry (2013) 3, e287; doi:1 0.1 038/tp.201 3.57; published online 23 July 2013 



Introduction 

Identifying genetic variants that increase the risk for a 
heritable disease may not only reveal the underlying biological 
mechanisms but might also provide new treatment targets. In 
one of the largest brain-imaging studies ever, a multinational 
consortium (ENIGMA, enhanced neuroimaging genetics 
through meta analysis) recently identified a common genetic 
variation in rs7294919 located on chromosome 12q24 being 
associated with hippocampal volume (n = 21,151; P=6.7 
x 10" 16 ): 1 ' 2 an increased number of rs729491 9-7" risk alleles 
was associated with decreased hippocampal volume. Hippo- 
campal volume is highly heritable and implicated in many 
(neuro)psychiatric disorders, such as schizophrenia, bipolar 
disorder, major depression and Alzheimer's disease. The 
hippocampus is a key brain structure for memory formation 
and stress regulation, brain functions that are strongly 
implicated in the above mentioned disorders. The identified 
quantitative trait locus for hippocampal volume may exert its 
effects by regulating the expression of the nearby TESC 
gene. 1 This previously not very well studied gene is expressed 
during development and also moderately during adulthood. Its 
protein product is involved in cell proliferation and differentia- 
tion and thus relevant for hippocampal integrity. 2 ' 3 



In the study reported here, we explored whether this 
quantitative trait locus also affects hippocampal function, 
using a previously established neuroimaging paradigm. 4 As 
the reliability of our target measure is an important issue in 
estimating specific genetic effects on hippocampal function, 
we further performed a test-retest reliability analysis of the 
identified contrasts to quantify the stability of our paradigm 
over time. 5 

Methods 

Subjects. A total of 224 German volunteers with grand- 
parents of European ancestry were recruited at Bonn, 
Mannheim and Berlin and genotyped for rs7294919. The 
study was approved by the local ethics committee and all 
participants gave prior written informed consent. No partici- 
pant reported a lifetime or family history of (neuro)psychiatric 
disease. Of our sample, 176 participants were rs7294919 TT 
homozygotes, 42 were CT heterozygotes, and 6 were CC 
homozygotes. No deviation from Hardy-Weinberg equili- 
brium was observed (x 2 = 0.027, df = 1 , P= 0.87). There was 
no genotype effect on demographic variables, behavioral or 
psychometric data (Table 1). 
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Table 1 Sample characteristics and behavioral measures 





TT 


CT 


CC 


P 


Subjects 


176 


42 


6 




Site 


79/56/41 


20/14/8 


5/0/1 


0.38 


(Bonn/Mannheim/ 










Berlin) 










Sex (male/female) 


81/95 


24/18 


1/5 


0.14 


Handedness 


159/12/3 


37/4/0 


5/1/0 


0.89 


(right/left/both) a 










Age (years) 


34.1 ±10.0 


34.8 ±9.9 


34.1 ±10.0 


0.51 


School education 


6.2 ±1.0 


6.2 ±1.3 


6.0 ±1.1 


0.87 


(years) 










Hits recall 


11. 6 ±1.8 


11. 4 ±1.9 


9.8 ±1.7 


0.07 


VLMT early recall 


12.4±2.4 


13.3±1.7 


13.0±2.5 


0.02 


VLMT late recall 


12.6±2.4 


13.8± 1.6 


13.7±1.4 


0.004 


VLMT recognition 


14.4±1.0 


14.6±0.7 


15.0±0.0 


0.19 


MWTB 


30.4 ±3.5 


30.9 ±2.8 


29.2 ±2.5 


0.46 


BDI 


2.8 ±3.3 


2.9±3.1 


2.4 ±2.9 


0.95 


STAI-S 


31 .3 ±0.6 


31 .6 ±1.2 


28.5 ±2.7 


0.95 


STAI-T 


35.2 ±0.8 


35.5 ±1.7 


36.4 ±3.9 


0.56 


Hippocampal 


4535 ± 473 


461 5 ±454 


4441 ±414 


0.53 


volume right (mm 3 ) 










Hippocampal 


4440 ±444 


4579 ± 426 


4466 ± 436 


0.19 


volume left (mm 3 ) 










Abbreviations: BDI, Beck Depression 
Mehrfachwahl-Wortschatz-lntelligenztest; 


Inventory; 
STAI-S, 


C, cytosine; 
State Trait 


MWTB, 
Anxiety 



Inventory — state version; STAI-T, State Trait Anxiety Inventory — trait version; 
T, tyronine (risk allele); VLMT, verbal learning and memory test. 
a Data missing for three subjects. 
Bold values indicate significant P<0.05. 

DNA extraction and genotyping. EDTA anti-coagulated 
venous blood samples were collected from all individuals. 
Lymphocyte DNA was isolated using the Chemagic Magnetic 
Separation Module I (Chemagen, Baesweiler, Germany). 
The single-nucleotide polymorphisms rs7294919 were geno- 
typed using a TaqMan 5' nuclease assay. 

Functional imaging task. During functional magnetic reso- 
nance imaging (fMRI), subjects completed an episodic recall 
task that was part of three consecutive memory tasks, that is, 
encoding, recall and recognition of face-profession pairs. 
These tasks were part of a functional imaging genetics 
battery. 6 The encoding task consisted of 16 face-profession 
pairs and 24 head contours as control condition with 4 blocks 
of 4 face-profession pairs and 4 blocks of 6 head contours 
each. Face-profession pairs were presented for 6 s, and 
head contours for 4 s. Each block lasted 24 s. Participants 
had to imagine the person acting in a scene of the written 
profession and had to indicate whether the profession suited 
the presented face. During the control condition participants 
had to indicate which ear of the depicted head contour was 
larger. The alternating sequence of four face-profession 
association blocks and four control blocks was presented 
twice in order to ensure successful encoding. During recall, 
the task in focus of the present report, faces were presented 
together with the question whether the depicted person had 
to complete apprenticeship or academic studies in order to 
qualify for the respective profession that had been learned 
during encoding. Subjects had to indicate by button press, 
which qualification was correct. Stimulus duration and control 
condition were similar to encoding, blocks were presented 
only once. For recognition testing faces were depicted 
together with two written professions and subjects had to 



indicate which profession was correct. Stimulus duration for 
recognition was 3 s. The control condition consists of 4 
blocks of 4 head contours each (3 s). Thus, each recognition 
block lasted 12 s. Similar to recall, blocks were only 
presented once. 

Behavioral and psychopathological measures. On a 

second day, subjects underwent neuropsychological assess- 
ment including testing for verbal intelligence (MWTB, 
Mehrfachwahl-Wortschatz-lntelligenztest) and memory 
(VLMT, Verbaler Lern und Merkfaehigkeits-Test, translated 
verbal version of Rey auditory learning task) and completed 
versions of the Beck Depression Inventory (BDI) and the 
State Trait Anxiety Inventory (STAI, Trait version). 

Imaging parameters. BOLD fMRI was performed on a 
Siemens Trio 3T scanner at the Life and Brain Center, 
University of Bonn, the Central Institute of Mental Health, 
Mannheim and the Charite University Hospital Berlin (para- 
meters: 33 slices, axially tilted (-30°), slice thickness 
2.4 + 0.6 mm gap, FOV 192 mm, TR 1.96 s, TE 30 ms, flip 
angle 80°). High resolution three-dimensional T1 -weighted 
images were acquired with 160 contiguous sagittal slices of 
1 mm thickness (FOV 256 mm, TR 1.57 s, TE 3.42 ms, Tl 
800 ms, flip angle 15°). 

Structural image processing. Structural MRI data of this 
sample were part of the ENIGMA sample. 1 The volumetric 
segmentation of the left and right hippocampus was 
performed with the Freesurfer image analysis suite (version 
5.0; http://surfer.nmr.mgh.harvard.edu/). The three mprage- 
images were fed into the volumetric stream which included 
motion correction, averaging, removal of non-brain tissue, 
automated Talairach transformation, intensity nomalization 
and segmentation of the subcortical white matter and deep 
gray matter volumetric structures (including the hippocam- 
pus) with the help of a probabilistic atlas. 7 ' 8 Quality control 
included image inspection, outlier detection via boxplot and 
subsequent manual editing of the brainmask, details are 
given in Stein et al. 201 2 2 

Functional image processing. Image processing and 
statistical analyses were conducted using statistical para- 
metric mapping methods as implemented in SPM8 (http:// 
www.fil.ion.ucl.ac.uk/spm/software/spm8/). Briefly, images 
were realigned to a mean image (movement parameters 
were confined to <3mm translation and <3° rotation 
between volumes), slice time corrected, spatially normalized 
to a standard stereotactic space (a brain template created by 
the Montreal Neurological Institute) with volume units 
(voxels) of 2 x 2 x 2 mm 3 , smoothed with an 8 mm full-width 
at half maximum Gaussian filter and ratio normalized to the 
whole-brain global mean. A first-level fixed-effects model 
was computed for each participant. Regressors were created 
from the time course of the two experimental conditions 
(memory, control) and convolved with a canonical hemody- 
namic response function. Movement parameters were 
included in the first-level model as regressors of no interest. 
For each subject, statistical contrast images of memory 
versus control were obtained. To test for genetic association, 
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these were analyzed using the general linear model in a 
second-level multiple regression analysis. To test for genetic 
association, these were analyzed using the general linear 
model in a second-level multiple regression random effects 
analysis with age, gender, scanner site and hippocampal 
volume as nuisance covariates to identify genotype effects 
on functional activation. 

Statistical inference. For all imaging methods, the signifi- 
cance threshold was set to P<0.05, corrected for multiple 
comparisons within our a priori defined anatomical region of 
interest, the hippocampus. For all analyses, we employed 
conservative analysis statistics by using family wise error 
rate, that exert strong control of type I error over multiple 
comparisons in imaging genetics. Regions of interest were 
defined a priori using anatomical labels provided by the 
Wake Forest University PickAtlas (www.fmri.wfubmc.edu/ 
downloads). 

Behavioral and neuropsychological data were analyzed 
using PASW Statistics 20 (IBM SPSS Statistics, Chicago, IL, 
USA) and x 2 -tests (site, sex, handedness), univariate (age, 
education, MWTB, performance, BDI) as well as multiva- 
riate (VLMT, STAI) ANOVAs with site, age and sex as 
covariates. 

Test-retest reliability of hippocampal activation. In order 
to quantify the robustness of our activation measures we 
analyzed test-retest reliability in an independent cohort of 
young healthy subjects (n = 20, 11 females, mean age 24.3, 
mean retest interval 14.5 days) using intraclass correlation 
coefficients (ICC) as defined by Shrout and Fleiss, 9 that is, 
the absolute (ICC(2,1)) and relative (ICC(3,1)) agreement of 
observations between sessions (for details see Plichta 
et al. 5 ). For the purpose of this study, the analysis was 
spatially restricted to the hippocampal cluster that exhibited 
significant genotype-dependent activation differences. Out- 
come measures of interest were ICCs of mean BOLD 
amplitudes of voxels within this mask. Analyses were 
performed using PASW Statistics 20 (IBM SPSS Statistics). 
Consistent with the prior literature, we interpreted ICC values 
below 0.4 as indicative of poor test-retest reliability, while 
ICC values between 0.4 and 0.75 were interpreted as 'fair to 
good' and those above 0.75 as 'excellent', respectively. 




right Hippocampus 

at [20, -16, -16] 



o 
o 




TT CT CC 

Figure 1 Carriers of the of the risk variant (T) for rs729491 9 exhibit significantly 
increased allele-dosage-dependent activation of the right hippocampus during recall 
(Z=3.45, P<0.05 family wise error (FWE) corrected for multiple testing across 
region of interest). Each red dot represents size of effect in one subject. 



Results 

We observed an effect of rs7294919 genotype in the right 
hippocampus (x=20, y= 16, z= 16, Z=3.45) with 
activation increasing with the number of T alleles (Figure 1). 
Further, in the external VLMT recall test performance 
decreased significantly with number of T alleles (see 
Table 1 ). Rs729491 9 had no impact on fMRI task performance, 
verbal intelligence, or depression and anxiety measures. 

Hippocampal volume did not differ significantly between 
genotype groups in our (sub)sample. This might be due to the 
smaller sample size compared with the sample used in the 
ENIGMA consortium. 

The test-retest reliability analysis provided evidence for 
above average reliability 10 of our recall task and a very good 
robustness of the hippocampal activation estimates that 



yielded significant differences in the present study (Intraclass 
correlations: ICC(2,1) = 0.73 (CI: 0.31 0.89), ICC(3,1) = 0.73 
(CI: 0.31 0.89)). 

Discussion 

Our findings implicate the variant identified by the ENIGMA 
consortium at rs7294919 in hippocampal function and 
episodic memory recall. This underlines the idea that 
rs729491 9 may affect the individual capacity to resist disease 
in terms of diminishing or boosting hippocampal resources. 
A test-retest analysis demonstrates a high robustness of our 
activation estimates in the hippocampus, both in relation to 
previously established formal criteria 9 and to the mean 
reliability of fMRI studies in general. 10 
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The observed increase of functional activation during 
memory recall in carriers of the rs27491929-7 risk variant 
might indicate a compensatory mechanism: as genotype 
groups did not differ in task performance during fMRI, 
increased hippocampal activation can be interpreted as 
overrecruitment needed for similar performance. 11 Reduced 
performance in an unrelated and more demanding VLMT 
might in turn indicate reduced reserve capacity when task 
difficulty is high. This diminished neuronal reserve could be 
related to the previously detected association between 
rs7294919 risk allele and decreased hippocampal volume. 2 

An alternative explanation for increased hippocampal 
activation could be that hippocampal hyperactivity corre- 
sponds to a relative decrease in synaptic inhibition 12 through 
neuronal excitotoxicity. 13 Overactivation of synaptic receptors 
by excitatory neurotransmitters, for example, glutamate, can 
lead to high levels of calcium influx into the cell that, if not 
sufficiently buffered, can activate a number of enzymes 
leading to structural cell damage. 14 Interestingly, the protein 
product of the TESC gene, tescalcin, is a calcium binding 
protein 15 and the expression of TESC has been associated 
with the rs7294919 genotype. 2 However, the exact nature of 
the relation between TESC expression and hippocampal 
volume and function remains to be determined in further 
translational studies and results of the present study need 
independent replication. 
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